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Abstract 

We review recently developed two-step approach for description of electronic and magnetic prop- 
erties of a new class of materials, the diluted magnetic semiconductors. In the first step we con- 
struct, on the basis of the state-of-the-art first-principles electronic structure calculations, the 
effective Ising and Heisenberg Hamiltonians which describe, respectively, the alloy phase stability 
and the magnetic excitations in the system. In the second step, we analyze properties of these 
effective Hamiltonians by various methods of statistical physics. As a case study, the prototypical 
diluted magnetic semiconductor Gai_ x Mn x As is studied in detail. We determine, among others, a 
possibility for clustering in this system, formation energies of various compensating defects, and es- 
timate short-range order parameters describing ordering tendencies in a system. On the other hand, 
by using recently developed local random-phase approximation approach, we evaluate the system 
Curie temperature and demonstrate its strong dependence on the sample preparation. We also 
emphasize the relevance of proper inclusion of the randomness in positions of magnetic impurities 
for a reliable estimate of the system critical temperature. Finally, we compare calculated Curie 
temperatures with available experimental data and briefly mention relation to other theoretical 
approaches. 
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I. INTRODUCTION 



Properties of semiconductors are very sensitive to a small amount of various impurities 
and defects while magnetism is a collective phenomenon often stable to high temperatures. 
Magnetic order in the system strongly influences other material properties, like e.g., the 
phase stability, and transport and optical properties. Also, unique properties of both semi- 
conductors and magnetic materials form basis of many important technologies. The com- 
bination of properties of both these classes of materials seems to be very promising and in 
fact resulted recently [l| into the discovery of a new class of materials, the diluted magnetic 
semiconductors (DMS), of which the diluted III-V DMS is the best example and (Ga,Mn)As 
system is the most frequently studied one both experimentally and theoretically. 

In addition to being promising in future spintronics applications, the DMS represent also 
the challenge to the solid state and material science physics. A main feature of these new 
materials is ferromagnetism which is primarily due to coupling of impurity magnetic mo- 
ments mediated by free carriers, typically by holes, in the host semiconductor valence band. 
These holes are introduced into the GaAs host by Mn-impurities but their concentration 
in the system can be strongly influenced by the presence of native defects, namely by As- 
antisites on the Ga-sublattice and by Mn-interstitials. The hybridization between magnetic 
impurity and host electronic states leads to a formation of virtual bound states for which 
proper inclusion of electron correlations can be quite important. However, one of the most 
important features characterizing the DMS is the presence of the disorder and, in particu- 
lar, how the disorder influences the distribution of magnetic and other impurities: the atom 
clustering or segregation, the tendency of foreign atoms to be incorporated into the system 
in the presence of other impurities, both native and doping ones, etc. are relevant questions 
to be addressed. In addition to these basic structural issues there is an equally important 
question of how important can be disorder for the Curie temperature of the system: suffi- 
ciently robust ferromagnetism is a key point for any future technological application of the 
DMS. 

There is an extensive literature on the theory of DMS which, however, summarizes 

mostly results of model approach to the problem. The parameter-free, first-principles studies 
are less frequent and a systematic, comprehensive review on this subject is still missing. The 
aim of this paper is not to give such a review but rather to illustrate one specific approach 
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which employs the first principles calculations as a starting point for the construction of 
simple, effective Hamiltonians which could address various important problems from the 
fields of the phase stability and magnetic excitations in a simple, transparent way rather 
than to attempt to study these problems fully on the first-principle level, which is still in 
most cases numerically prohibitive. Below we will consider two such models, the Ising and 
Heisenberg Hamiltonians, whose parameters, the effective pair interactions in the former 
case, and the pair exchange interactions in latter case, are determined on the basis of the 
same first principles calculations. In both cases such parameters can be determined either 
from ad hoc chosen structural or magnetic configurations in the framework of conventional, 
electronic structure supercell calculations 0, Q or, as in the present case, by explicitly 
including the effect of disorder in terms of the coherent-potential approximation (CPA) |7|. 
One obvious advantage of the present approach is the possibility to include the effect of small 
concentrations of various types of impurities as well as the effect of finite lifetime due to 
disorder, features typical for the DMS, on the same footing. On the other hand, the effect of 
clustering on the electronic and magnetic structure can be more straightforwardly included 
using the supercell approach. A detailed comparison of both approaches is, however, beyond 
the scope of the present paper. An emphasis will be put on the study of magnetic properties 
of GaMnAs alloys, in particular on the determination of their Curie temperatures based on 
a detailed study of the corresponding Heisenberg Hamiltonian by sophisticated statistical 
methods which take into account both random positions of Mn-impurities and the presence 
of compensating defects in a system. 

II. ELECTRONIC STRUCTURE 

We have determined the electronic structure of the DMS in the framework of the first 
principles all-electron tight-binding linear muffin-tin orbital (TB-LMTO) method in the 
atomic-sphere approximation using empty spheres in interstitial tetrahedral positions of the 
zinc-blende lattice for a good space filling. We used equal Wigner-Seitz radii for all atoms 
and empty spheres. The valence basis consists of s-, p-, and <i-orbitals, we include scalar- 
relativistic corrections but neglect the spin-orbit effects. The substitutional disorder due to 
Mn-atoms and other possible defects is included within the CPA. The charge selfconsistency 
is treated in the framework of the local spin density approximation using Vosko-Wilk-Nusair 
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parameterization for the exchange-correlation potential 8|. The lattice constant of the pure 
GaAs (a = 5.653 A) was used in all calculations but we have verified that we can neglect 
a weak dependence of the sample volume on defect concentrations. Further details of the 
method can be found in Q]. 

The magnetic disorder is treated in the framework of the disordered local moment (DLM) 
method 9j which is the simplest way of including disorder in spin orientations and which 
is justified for atoms with large exchange splitting. A comparison of the total energies of 
the ferromagnetic (FM) state with its DLM counterpart is the simplest way to investigate 



the magnetic stability of the DMS alloy 10, • In addition, the DLM is a natural starting 



point for the phase stability studies due to the fact that sample preparation is done at 
temperatures well above the system Curie temperature. The DLM can be included in the 
framework of the CPA: the Mn atoms have collinear but random positive (Mn + ) and negative 
(Mn~) orientations with corresponding concentrations x + and x~, x = x + + x~ , where x 
is the total Mn-concentration. The degree of magnetic order can then be characterized by 
the order parameter r = (x + — x~~)/x, and x ± = (1 ± r)x/2. In the (FM) state, r = 1, all 
magnetic moments are aligned in the direction of a global magnetization. The non-magnetic 
state, r = 0, is characterized by a complete disorder of spin directions with vanishing total 
magnetization while a partial ferromagnetic state is characterized by < r < 1. For more 
details concerning the determination of the magnetic phase diagram and magnetic moments 
in different magnetic states we refer to fl^l . [Lsl ] . 

In Fig. [H we plot the local density of states on Mn-atoms, clearly illustrating general 
features of the GaMnAs alloys which will be relevant for their magnetic properties: (i) the 
halfmetallic character of spin-subbands with a gap in the minority states, and (ii) the strong 
disorder introduced by the presence of Mn-impurities in the GaAs host as manifested by 
its pronounced difference from the host GaAs density of states. Also the reference host 
bandstructure is shown for an illustration. 



III. ISING HAMILTONIAN OF AN ALLOY 

Some important problems concerning the structure and phase stability of DMS can be 
studied on the basis of the Ising Hamiltonian describing various configurations of disordered 
alloys with parameters determined from first principles. A particular configuration of a 
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FIG. 1: Electronic structure of (Gao.o5,Mno.o5)As random alloy. The upper and lower insets show 
the total density of states and the bandstructure along principal axis in the Brillouin zone of the 
host GaAs crystal. The spin-polarized local density of states on Mn-atoms substituting cations on 
Ga-sublattice is shown in the main frame. 

homogeneous disordered multicomponent alloy is characterized by occupation indices i]^, 
where 77^ = 1 if the site R is occupied by an atom of type Q, and 77^ = otherwise. 
Configurational averaging of occupation indices (?7^) = c Q yields the concentrations c Q . 
The configurationally dependent part of the alloy internal energy is given by the effective 
Ising Hamiltonian 

H = +\Y.Y. V ^V Q ^% + ..., (1) 

RR' QQ' 

where V^^, are interatomic pair interactions. Here we will not consider higher order inter- 
actions. 

Of)' 

The pair interactions V-^j in DMSs consist of two contributions, 

V QQ' _ ,.QQ' , aQQ' (9 \ 

RR' — RR' ~+~ rRR' j \ z ) 
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where the v^^, result from mapping of the band part of the total energy onto the Ising 
Hamiltonian ([1]) and the represent the electrostatic interaction energy of a pair of atoms 
Q, Q' located at sites R, R' (for derivation see Appendix A in [l4]) = e 2 g^g^/|R — R'|, 
where q® a = — q is the effective net charge of atomic species Q defined as a difference of the 



net charge <p of atomic species Q and the averaged charge q. The 
is calculated using the Generalized Perturbation Method (GPM)|5 



3and term contribution 
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,,QQ' 

RR' 



— Im 

7T 



dE tr 



(3) 



where tr denotes a trace over angular momentum indices (£m) and the spin index a, z = 
E+iO, E-p is the CPA Fermi energy, E min is a suitably chosen energy below the valence energy 
spectrum, <?rr' (z) denotes the block of the averaged auxiliary Green function between sites 
R and R', and t^(z) is the t-matrix for atomic species Q. 



A. Impurity formation energies 

The calculated total energies (per elementary cell) make possible to investigate impu- 
rity formation energies for DMSs. The formation energy e[Ab] of an impurity Ab which 
substitutes a host atom B in a (generally multicomponent) alloy A i .B 1 _ x is defined as 

e[A B ] = NE[A x+ g x Bi^ x -5 X ] + E at [B) - {NE[A X B^ X ) + E at [A]} , (4) 

where N is the number of elementary cells in the alloy, Sx — 1/N, -Eat [A] is the energy of 
an isolated atom A, and EfA^Bi-j;] is the energy of the alloy per one elementary cell. By 
expanding into linear terms in 5x one finds 

e[A B ] = ggMlzd + Eat [B] - E«[A] . (5) 

Our calculations [2| for (Gai_ I __j / Mn ;! ;Asj / )As alloys have shown that the formation energy 
of As antisite defect decreases with increasing concentration x of substitutional Mn atoms, 
which means that the number of the antisite defects can be considerably enhanced in the 
presence of substitutional Mn. Similarly, the formation energy of the substitutional Mn 
decreases with an increasing concentration y of As antisites. This means that the presence 



of As antisites (and probably also of other donors 16[]) is important for an improved solubility 
of Mn in III-V materials. 
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The behavior of formation energies for interstitial Mn impurities and for As antisites in 
the presence of interstitials is opposite to that found for substitutional Mn because increasing 
concentration of one of the species leads to a growth of the impurity formation energy of 
the other one. This is in agreement with the growth mechanism of (Gai_ rr _ y Mn a ,As 2y )As 
alloys discussed by Erwin and Petukhov [17|: the Mn atoms are first incorporated into 
interstitial positions under low concentration of As antisites, and later on, during growth 
or annealing, as substitutional impurities. In particular, the conversion of the interstitial 
Mn into a substitutional form is facilitated during the growth process if a sufficient number 
of As antisites is available. Both discussed mechanisms contribute to the self-compensation 
behavior of as-grown (Ga,Mn)As alloys. 



B. Alloy stability 

The calculated total energies (per elementary cell) make possible to investigate the sta- 
bility of DMSs with respect to segregation into systems with different chemical composition. 
For example, consider segregation of (Gai_ x _ ?; Mn a .As y )As into an alloy without As-antisites 
(Gai_ x Mna;)As and an alloy with the highest possible concentration of As-antisites which is 
still not overcompensated, (Ga!_ 3;c /2Mn x As x /2)As, 

AE(x,y) = ^(Ga^^Mn.As^As] - ^-^[(Ga^Mn^As] 
2v 

-JLE[{G^ 3x/2 Mn x As x/2 )As} . (6) 



The energy AE(x, y) is negative 14( which indicates the stabilizing effect of the As antisites. 
This is in agreement with the above conclusions based on the impurity formation energies. 

Similarly, we can consider segregation of the alloy (Gai_ I Mn a ,)As into pure GaAs and an 
alloy with higher concentration of Mn atoms, say xo, (Ga.i- XQ Mia Xo )As, 

AE(x,x ) = E[(Ga 1 _ a Mn g )As] - X ° ~ X ff [GaAs] 

_^ jB [(Ga 1 _ a;o Mn xo )As] . (7) 

X Q 

Our calculations have shown that this energy is positive which shows that the alloy 
(Ga 1 ^ a .Mn 3 .)As is thermodynamically unstable with respect to segregation. 
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C. Ordering tendencies 



The effective interatomic pair interactions in semiconductors decrease rather slowly 
with the interatomic distance |R — R'| and have to be calculated over many coordination 
spheres 141 ] . Also their Coulombic part, even though weak, is long-ranged. In order to 
analyze possible ordering patterns, we employed the linearized concentration wave method 
and for temperatures above the ordering temperature T or( j we also calculated Warren-Cowley 
short-range order parameters. One has to keep in mind that in these methods only the con- 
figurational part of the entropy is taken into account which usually leads to an overestimation 
of ordering temperatures. 

The ordering temperature T ord and the type of ordered structure that appears below 
T or d can be studied in terms of the concentration- wave method 5), [lsj. Here we employ 
its linearized version [3] extended to a multicomponent alloy since we consider possible 
ordering of four atomic species (Ga, Mn^, Mn^, As) on the cation sublattice. In a mean field 
approximation (i.e., assuming a Bragg- Williams form of the entropy {20)]) the free energy is 
expressed in terms of local concentrations c^, 

F = \ E E + * bT EE . (8) 

RR' QQ' R Q 

where k-g is the Boltzmann constant and T is the temperature. Starting from the disordered 

state, the free energy can be expanded up to quadratic terms in concentration fluctuations 
x c Q _ C Q _ r Q 



F = F ° + ^EE [ V nS + k ~f^S QQ ] 64 5c% (9) 



1 Sy^ - ^ V ^ i", -t )( )' /•'!!•/ 
RR' 

with terms linear in <5c^ vanishing because Y^, = for all Q, and Ylwa 1 ^rr' c<3 ' is a 
constant for all R and Q (for details see Ref. llSj). Here Fq is the free energy in the absence 
of concentration waves (c^ = for all R). Equation (Q can be rewritten in terms of a 
lattice Fourier transform in a matrix notation as 

BZ 

AF = F - F = - ^(k) [V(k) + k*TC- 1 ] E(k)|e(k)| 2 , (10) 

k 

where [V(k)] QQ/ = V QQ '(k), and the matrix C is defined as [C] QQ > = c Q 5 QQ >. In Eq. (JTDP 
the concentration fluctuations <5c^(k) are expressed in terms of a vector E(k) and the order 
parameter e(k) as <5c^(k) = E^(k)e(k). At sufficiently high temperatures AF is positive 



definite, because the hermitian matrix Vity + k^TC^ 1 has only positive eigenvalues and thus 
the high temperature state is completely disordered (e(k) = for all k). With decreasing 
temperature it can become indefinite at T orc j because of a vanishing eigenvalue for a critical 
vector k which determines the period of the concentration wave. The components Y®(k) 
of the critical eigenvector determine the amplitude of the concentration wave for each alloy 
component Q. For each k, the minimization of AF, and thus the eigenvalue problem, is 
subject to the subsidiary condition X1q^ Q (^) = which follows from ^2qSc^ = 0, valid 
for each R. The ordering temperature is then found as the largest eigenvalue of the matrix 

0(k) = -k^cW^c*. 

For example, for a ferromagnetic alloy (Gao.93MnJ 06 Aso.oi) As we found [14] the ordering 
temperature T or d = 775 K and the ordering vector k = 0.274(1, 1, 1) /a, where a is the 
lattice constant. A closer examination shows that the largest eigenvalues of the matrix 
@(k) have very similar values (within 1 K) for k- vectors close to a surface of a sphere of 
radius 0.479/a, which corresponds to a domain of a characteristic radius 3.7 nm in real 
space. The components of the eigenvector Y(k ) are F(Ga) = 0.726, F(IvW) = —0.686, 
and K(As) = —0.040, which corresponds to formation of domains of two types: in the first 
type the concentration of impurities Mn and As is increased, while in the second type the 
impurity concentrations are diminished in agreement with the above estimate based on the 
segregation energy ([7j). 



D. Short-range order parameters 

The Warren- Cowley short-range order parameters 

ER '~~ {>&><>&) <f><f 1 ' 

provide detailed information on mutual correlations of impurities; they also can directly be 
used in calculations of transport and magnetic properties. The matrix of the Warren-Cowley 
parameters can be approximately calculated by means of the Krivoglaz-Clapp-Moss (KCM) 
formula in reciprocal space, 

a(k) = —D[M + /#(k)] _1 D T , (12) 

where (3 = (fc B T) _1 , and the matrix M is defined as [M]qq/ = l/c° + <5q,q' /c q '. The matrix 
D is introduced to ensure a correct normalization of a for R = R' which follows from the 
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definition in <^2>. Note that V(k) is a Fourier transform of = +V^ R ,-V$! L ,-V£% 
, which is the matrix of pair interactions from which one atomic species was eliminated. The 
inverse lattice Fourier transform of (fl2l) then yields the Warren-Cowley parameters in the 
real space. 

The results for the Warren-Cowley parameters show a strong tendency to an aggregation 
of Mn atoms with the same orientation of spin. This finding is in line with the results of Ref. 
2l| . Close pairs of Mn atoms with opposite moments become possible, while the probability 
to find close pairs of As antisites is very low. A moderate aggregation of Mn atoms and As 
antisites is also possible. 

The results for the Warren-Cowley parameters show a strong tendency to an aggregation 
of Mn atoms with the same orientation of spin. This finding is in line with the results of Ref. 



2l| . Close pairs of Mn atoms with opposite moments become possible, while the probability 



to find close pairs of As antisites is very low. A moderate aggregation of Mn atoms and As 
antisites is also possible. 



E. Conclusions from structural studies 



The main conclusions can be summarized as follows: 

(i) The alloys are thermodynamically unstable with respect to segregation into related 
compounds or alloys. 

(ii) As-antisites have a stabilizing effect and make the incorporation of substitutional 
Mn atoms energetically favorable. On the other hand, incorporation of Mn atoms into 
interstitial positions is energetically favorable only at low concentration of As antisites. 

(iii) Formation of domains of two types, namely, with an enhanced and with a lowered 
concentration of impurities (substitutional Mn atoms and As antisites), can be expected. 
The characteristic size of the domains depends on chemical composition and might be of 
order of several nm. 

(iv) A strong tendency to aggregation of substitutional Mn atoms with parallel magnetic 
moments is expected, while the formation of close pairs of As antisites is highly unlikely. 
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IV. HEISENBERG HAMILTONIAN 



The knowledge of exchange interactions allows one to address in detail the character of 
magnetic excitations in the DMS, i.e., to evaluate the Curie temperature, the spin-wave 
stiffness, and the spectrum of low-lying magnetic excitations. Magnetic excitations in ferro- 
magnets are of two different kinds, namely, Stoner excitations associated with longitudinal 
fluctuations of the magnetization, and spin-waves, or magnons corresponding to collective 
transverse fluctuations of the magnetization direction. The low-temperature regime is dom- 
inated by magnons and Stoner excitations can be usually neglected. In this section we 
will construct an effective random Heisenberg Hamiltonian with classical spins which will 
be used in the next section to study the critical temperatures. The main idea is to sep- 
arate a relevant part of (very) small magnetic energies responsible for magnetic ordering 



221. The 



from huge total energies obtained from first-principles total energy calculations 
validity of this approach, based on the adiabatic approximation, is in particular justified 
for magnetic atoms with large exchange splitting, like, e.g., Mn-impurities in GaAs host. 



The mapping is further simplified by using the magnetic force-theorem [22j, [23| which states 
that the band energy of the calculated ground state of a reference spin structure can be 
used as an estimate for corresponding total-energy differences in the excited state. Note 
that intracell non-collinearity of the spin polarization is neglected since in this approach 
we are primarily interested in low-energy excitations due to intercell non-collinearity. The 
application of this approach to disordered systems like the DMS is significantly simplified by 
using the vertex-cancellation theorem (VCT) which justifies the neglect of disorder- induced 
vertex corrections in Eq. (THl) below. The VCT was derived in 24| under rather general con- 
ditions and facilitates an efficient evaluation of exchange interactions, exchange stiffnesses, 
spin-wave energies, etc. The present real-space approach is particularly suitable for random 
systems with low concentrations of magnetic impurities such as the DMS, where the effect 
of disorder is treated in the framework of the CPA. It is also possible to estimate exchange 
interactions from energy differences between parallel and various antiparallel configurations 
of a few magnetic atoms in a supercell representing specific concentration (the real-space 
supercell approach) [25|. In this way exchange interactions of magnetic clusters could be 
estimated approximately. Alternatively, the reciprocal-space approach to the problem is 
also possible, namely the frozen-magnon method combined with the supercell approach 26]. 



II 



We refer reader to a recent review [6| for more details, comparison of the real-space and 
reciprocal-space approaches, as well as for applications of the above formalism to a broad 
range of various magnetic systems. 



A. Effective pair exchange interactions 

The mapping of the total energy of itinerant electron system connected with small rigid 
rotations of two magnetic moments at sites R and R' can be described by the effective 
Heisenberg Hamiltonian with classical spins as discussed above 

H eS - ^ JR,R/e R ■ e R / . (13) 

R^R 

Here, Jr,r' is the exchange interaction energy between sites (R, R'), and e R , e R are unit 
vectors pointing in directions of local magnetic moments at sites (R, R'), respectively. In 
the present formulation the values and signs of the magnetic moments are already absorbed 
in the definition of the J^, so that posifve (negative) J^, correspond tofe^ag- 



netic (antiferromagnetic) coupling. By adopting the magnetic force-theorem y, |22|, |23j , the 



configurationally averaged effective pair exchange interactions J RR/ between two magnetic 



atoms M, M' located randomly at sites R and R' are given by the following expression [ll| : 

J r,r = ^ Im J tT L [tin ( z ) ^r,r ( z ) V ( z ) 9 R >; R dz . (14) 

Here, tr^ denotes the trace over angular momenta L = (£m), the energy integration is 
performed in the upper half of the complex energy plane along a contour C starting below the 
bottom of the valence band and ending at the Fermi energy, and 8^(z) = P^^(z) — P^^(z), 
where the P R ' a (z) are the L-diagonal matrices of potential functions of the TB-LMTO 
method for a =|, \ corresponding to a particular magnetic atom M. The matrix 8^(z) 
reflects the exchange splitting of atom M. The quantities ^ RR f \z) and g^, R \z) refer to 
site off-diagonal blocks of the conditionally averaged Green function [7] , namely, the average 
of the Green function over all configurations with atoms of the types M and M' fixed at 
sites R and R', respectively, determined in the framework of the CPA [7j. The exchange 
interactions between magnetic moments induced on non-magnetic atoms are negligible as 
compared to the exchange interactions between magnetic atoms in the present case, i.e., 
M=M'=Mn in (Ga,Mn)As alloys. The main advantage of the present approach is the 
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explicit expression for J^^f which can be evaluated straightforwardly even for large distances 
d = |R— R| between sites R and R' and thus allowing the study of their asymptotic behavior 
as a function of the interatomic distance d [s, It should be noted that the rigidity of 
spins during rotations as required by the Heisenberg model is preserved by construction 
in the present approach and in the reciprocal-space method [26| but not in the real-space 
supercell approach [25]. 

The effect of impurities on the host bandstructure is usually neglected in model theories 



3|, but it need not be a small perturbation in the presence of the virtua 
271 ] . It is also relevant for DMS systems as it was demonstrated recently 28 



-bound states 
by using the 



selfconsistent local RPA theory 29|: the combined effect of random geometry and thermal 
fluctuations is crucial and calls for a proper treatment. In particular, it was shown that 
damped RKKY interactions often used in model studies, are unable to represent properly 
the effect of virtual bound states in the host band on the values of exchange interactions, and 
corrections are needed (see also 30j). The standard RKKY exchange interactions cannot 



explain ferro magnetism observed in the DMS: one has to go further in this perturbation 
scheme to include properly effects of resonances due to virtual bound states. On the other 
hand as already mentioned, in the present approach is the effect of virtual bound states 
included by construction. 

The neglect of the effect of impurities on exchange interactions means that the unper- 
turbed host Green function appears in ()14p rather than its conditionally averaged counter- 
part (jftw a (z). The neglect of the renormalization of the host Green function by scatterings 
on impurities is two-fold: it introduces a phase factor and modifies the amplitude of the os- 
cillations as compared to the conventional RKKY formula 31]. 

The present approach neglects the effect of local environment effects: while the individual 
Jr fl' can ^ e very different for a particular environment, the corresponding configurationally 
averaged J^n' * n ran dom systems are close to the CPA value, Eq. (fill) , as demonstrated 
recently 32J . The CPA also correctly describes the concentration trends and the carrier con- 
centrations which, in turn, determine the size of the alloy Fermi surface and thus the periods 
of oscillations. The exchange interactions in GaMnAs alloys are exponentially damped due 
to the alloy disorder and by their halfmetallic character as illustrated in Fig. [2] below. The 
exchange interactions are also strongly anisotropic in the real space due to the underlying 
zinc-blende lattice while an additional, but weaker anisotropy, could be due to the neglected 



13 



2.5 
2 

1.5 
1 

0.5 


-0.5 



0.5 



2.5 



-5 
-7.5 
-10 









Mn 05 As y )As 


»- o- 








" . y=0.0 


o y=0.01 


• - 

V - 



m 



123456789 10 11 
(d/a) along [1 10] -direction 



1.5 2 
(d/a) 



2.5 



3.5 



FIG. 2: Exchange interactions j Mn > Mn between pairs of Mn atoms in (Gao.95- y ,Mno.o5, As J/ )As 
alloy plotted as a function of their interatomic distance d (in units of the lattice constant a). Full 
and empty symbols correspond to y = and y = 0.01, respectively, where y is the concentration 
of As-antisites on the Ga-sublattice. In the inset we show ln\(d/a) 3 J Mn ' Mn (d)| as a function of the 
interatomic distance d along the [110]-direction with dominating values of exchange interactions. 

spin-orbit effects. Finally, because the hybridization between Mn- and host-atoms is in- 
cluded by construction to all orders, both the RKKY-like and superexchange interactions 
are included in the theory although their separation is not possible in a simple way. 

The typical results for exchange interactions in (Ga,Mn)As alloys are illustrated in Fig. [21 
We observe well pronounced ferromagnetic behavior of GaAs semiconductor doped with 
5% of Mn- impurities: exchange interactions between pairs of Mn atoms, j Mn > Mn ; are fer- 
romagnetic over number of nearest neighbors with the dominating first nearest-neighbor 
interactions jj^ 11 ^ 11 but without standard RKKY-oscillations. Interactions are strongly re- 
duced by the presence of native As-antisite defects, which reduce the number of free carriers: 
while each Mn-impurity introduces one hole into the host valence band, the As-antisite (and 
similarly, the Mn-interstitial) introduces two electrons. With increasing concentration of 
As-antisites the leading j^ n ' Mn becomes negative (antiferromagnetic) and such frustration 
coincides with the extinction of the ferromagnetism in the system (see below). The damp- 
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ing of exchange interactions due to alloy disorder and, most importantly, due to halfmetal- 
lic behavior of GaMnAs alloy, is illustrated in the inset of Fig. [2j the linear decrease of 
amplitudes of the logarithm of exchange interactions multiplied by the RKKY-like factor 
(ln|(<i/a) 3 J Mn ' Mn |), indicating their exponential damping for large interatomic distances, is 
obvious. 



B. Curie temperatures of diluted (Ga,Mn)As alloys 

Here we briefly describe new approach for a quantitative determination of the Curie tem- 
perature T c which correctly takes into account the randomness in positions of magnetic 
impurities as well the presence of native defects, like As-antisites and Mn-interstitials and 
thus allows to explain their role in real samples. The novel feature as compared to the 
conventional magnet thermodynamics is the dilution, i.e., the presence of randomly dis- 
tributed magnetic defects of low, but finite concentration leading to the effect of magnetic 
percolation. It was clearly demonstrated recently JsSj that a sophisticated treatment of spin 
fluctuations using the RPA or Monte-Carlo methods but without inclusion of disorder is 
itself unable to explain experimentally observed Curie temperatures in the framework of the 
parameter-free approach: calculated critical temperatures are too high as compared to the 
experiment. Recently, few groups have formulated parameter-free approaches which take 
into account properly both disorder and jspin-fluctuations and arrived at a fair agreement 



with experiment for both well-annealed 29 



34 



35 



36, 



371 ] and as-grown samples [391]. 



The present accurate semi-analytical method separates the exact treatment of disorder 
and the RPA treatment of spin-fluctuations and has three stages: (i) the first-principle 
determination of exchange parameters of random classical Heisenberg model [ll]; (ii) gen- 
eration of a sequence of random configurations on a disordered lattice (fee Ga-sublattice 
for (Ga,Mn)As alloys) by MC sampling technique; and (iii) for each configuration the ran- 
dom Heisenberg model is treated analytically within the RPA. As the lattice is random, the 
equations have to be solved numerically. 

In the third step, the Green function Gr.r' for impurity spins at sites R and R' satisfies 



(E-h e i)G RiW (E) = 2(e z n )8 RjW -(e z n )Y,Jn,K."Gn»,B.>(E) 



R" 



hf = ^J R ,R"(e R ). 



(15) 
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The quantity /iff is the local effective field acting on the spin at the site R, and (e R ) is the 
local magnetic moment at the site R normalized with respect to the magnetization averaged 
over all impurities. For a given temperature, G-r..h/{E) for impurity spins are determined 
following the procedure similar to that of Callen 38] resulting in the selfconsistent solution 



of the following set of equations for the Curie temperature (for more details see paper 39]) 



k B T c = 
An,n(E) 



Gr,r(£) 



3-/Vi m p ^ F R 



Ak, r (E) 
E 



dE. 



2tt 



-Im- 



A 



R 



A R = lim T ^ Tc (e R ) / e z 



(16) 



where A^ mp is the number of impurity sites and e z denotes the averaged value of e R over 
all impurity sites. At the end, the average over typically few hundred configurations of 
the system, each of them including a few hundredth of impurity sites were enough to give 
robust values of the Curie temperature (with a numerical accuracy of about few K). On the 
contrary, the MFA value as evaluated of the average lattice which neglects the randomness 
in positions of magnetic ions is simply /c B T c MFA = (2x/3) J^ R J ,r, where x denotes the 



concentration of magnetic-atoms 



Hi 



The calculated T c according Eqs. (fl~6l) are shown in Fig. [3] for the case of fully-annealed 
(Ga,Mn)As samples without native defects (uncompensated samples) and with essentially 
random distribution of Mn-atoms 29j]. In the same Figure we also show experimental 
results of Edmonds et al. [40], Matsukura 41 1. and Chiba [42]. The agreement between 
theory and experiment is very good, except for the single highest concentration (9 %). Our 
calculations thus suggests that for this specific concentration the annealing is complete. The 



M 



371 ] based on Monte- Carlo 



present results are in a good agreement with related studies 
simulations used for the treatment of both spin-fluctuations and disorder: the use of the 
same exchange parameters gives for 5 % 137 K as compared to the value of 125 K in the 
present approach. Related ap pro ach which is based on the Korringa-Kohn-Rostoker (KKR) 
Green functions developed in 35J gives 103 K. The present theory also correctly predicts an 
expected threshold (about 1.5 %) below which there is no ferromagnetism. It should be noted 
that the theoretical threshold for the occurrence of magnetism on the fee lattice and the first 



3,0. 



nearest-neighbor Heisenberg model (19 %) was also confirmed by numerical studies [31 
The calculated realistic exchange interactions, however, extends over several neighbor shells 
which reduces the effect of magnetic percolation in realistic alloys. On the other hand for 
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FIG. 3: Calculated Curie temperatures for (Gai_ x Mn x )As alloys as a function of Mn- 
concentration. Theoretical results are compared with available experimental data for both well 
annealed and as- grown samples. 

as-grown samples Q, Q] the Curie temperature is significantly reduced by native defects. 

An immediate question is what is the proportion of two such defects, Mn-interstitials and 
As-antisites in as-grown samples. In order to answer this question we have first investigated 
the effect of As-antisites. The results are presented in Fig. HI clearly showing a weak de- 
pendence of Curie temperatures on the effective carrier concentration 7. First, we observe 
ferromagnetic instabilities for smaller values of 7. The reason for such instabilities is that for 
decreasing carrier concentrations are exchange interactions increasingly dominated by the 
antiferromagnetic (super exchange) contribution leading to a frustration. The experiment 
(40I ] for as-grown samples with a nominal Mn-concentration of x = 0.067 shows, however, 
no ferromagnetic instability. Second, relatively very weak dependence of calculated T c on 7 
contradicts experiment which, on the contrary, shows a pronounced dependence on the car- 
rier concentration Q . These fact s seem to exclude As-antisites as a dominating mechanism 
for compensation in as-grown samples. 

We will now demonstrate that observed Curie temperatures of unannealed or compen- 
sated samples can be explained assuming that interstitials defects dominate [39]. Such a 
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FIG. 4: The dependence of calculated Curie temperatures for (Gai- x - y Mn x As y )As alloys for 
three different Mn-concentrations and for varying concentrations y of As-antisites on Ga-sublattice 
plotted as a function of the carrier density. The parameter 7 is the ratio of the carrier concentration 
n/j = x — 2y and Mn-concentration x. 



dominance agrees with experimental observation of Wolos et al. 43) and Wang 44J. We 
emphasize that a clear proof of the role of interstitials is still necessary as other techniques 
ike transmission electron micrography [45] or infrared absorption and positron annihilation 



46( seem to indicate a much higher concentration of antisites. On the other hand, recent 



theoretical study of Wu 



471 ] also seems to support the model of as-grown alloys described 



below. The relevant quantity is the compensation parameter 7 = n/,/iM n , where rih is the 
carrier density and % n is the nominal concentration of Mn-atoms. It should be noted that 
both 7 and xun are available from the experiment 401 ] . 



Recent first-principle calculations 



481 ] and Rutherford backscattering experiments 49] 



indicate that Mn-interstitials (Mnj) are preferably attracted Mn-substitutional atoms on 
Ga-sublattice (Mnca) forming pairs of spins with a strongly antiferromagnetic coupling. 
We assume that Mn/ are not completely random but form, with a high accuracy, bound 
singlet pairs whose effect on magnetically active ions is very small. The remaining active 
Mn-atoms with effective concentration x e s = £Mn(Ga) — 2xMn(i) which are not directly cou- 
pled to interstitials are assumed to be distributed randomly and interact via the effective 
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FIG. 5: Curie temperatures of GaMnAs alloys as a function of Mn-concentration. Note that exper- 
imental data (diamonds) are plotted for effective Mn-concentration of magnetically active atoms 
and corresponding effective 7 as found in experiment. Squares and triangles are Curie temperatures 
corresponding to Mn-concentrations 0.035 and 0.05, respectively, calculated for densities of holes 
that correspond to the experiment. The small circles (dashed line) correspond to uncompensated 
samples 7 = 1. 

exchange coupling corresponding to a measured carrier density (or, alternatively via the 
experimentally determined parameter 7 e g). The main features of the above model are also 
supported by recent first-principles calculations {47J]. Technically, one needs to determine 
exchange parameters, Eq. (TT4l) . but for concentration x and the effective number of carriers 
nh treated as independent parameters. This is not a straightforward task in the framework 
of the first-principle theory without adjustable parameters, but it could be achieved with 
help of co-doping by impurities which have negligible influence on electronic properties of 



the alloy at the Fermi energy which is relevant for coupling between impurities [11J. For 
example, the Zn-doping on Ga-sublattice doping can increase the number of carriers. Sim- 
ilarly, the doping by Se-atoms on As-sublattice or by As-antisites decrease the number of 
carriers. It should be noted that co-doping is used here as a purely calculational device to 
control carrier density while keeping calculations fully self consistent. 
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The results are summarized in Fig. [5] where, for comparison, also results for 7 = 1 
(see also Fig. [3]) are shown (small differences are due to an improved statistics). In the 
same figure are also given experimental data for nominal Mn-concentration of 0.067 % but 



for a different annealing corresponding to different effective Mn-concentrations [39J. We 
observe, that well-annealed samples of highest T c are in very good agreement with the 
calculated values for uncompensated samples (7 = 1). We remark that 7 = 1 curve can 
be accurately parameterized up to XMn = 0.1 by the curve T c w A(xmu — x c ) 1 / 2 , where 
x c = 0.0088 and A = 649 K. Deviation from 7 = 1 curve is small for intermediate T c but 
it becomes increasingly visible for as-grown samples. For example, T c of as-grown sample 
which corresponds to 7 e g = 0.52 and x e s = 0.035 agrees well with the calculated value 
(square symbol) . A similarly good agreement is obtained also for other, differently annealed 
samples. One can thus conclude that the present parameter-free theory is able to account 
for both well-annealed and as-grown (Ga,Mn)As samples. Based on a good agreement with 
experiment one can also conclude that the dominating mechanism for the reduction of the 
Curie temperature of as-grown samples as compared to well-annealed ones is the presence 
of Mn-interstitials. 



V. CONCLUSIONS AND PROSPECTS 



We have presented a unified approach to describe electronic, structural, and magnetic 
properties of diluted magnetic semiconductors based on the GaAs host. Our approach is 
based on the first-principles electronic structure calculations which take into account low 
concentrations of various defects present in the system as well as the finite lifetime due 
to such defects. The calculated total energies are then used to construct simple effective 
Hamiltonians, namely the Ising model and the (classical) Heisenberg model, with parameters 
which are obtained from parameter-free first principles calculations. 

We have shown that a combination of first principles electronic structure calculations 
with a relatively simple methods of statistical mechanics applied to model Hamiltonians can 
give a coherent picture of the phase stability and possible ordering as well as of the critical 
temperatures in diluted magnetic III-V semiconductors prepared under different conditions. 
We have clearly demonstrated that inclusion of the randomness in positions of Mn-impurities 
is relevant for a good quantitative agreement between theory and experiment for both as- 
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grown and well-annealed samples. Based on a good agreement with experiment one can 
conclude that the dominating mechanism for the reduction of the Curie temperature of 
as-grown samples as compared to well-annealed ones is the presence of Mn-interstitials. 

One has to be aware that the samples studied in experiment need not be in thermody- 
namical equilibrium, but rather are in a metastable state corresponding to a local minimum 
of the thermodynamical potential. Phase stability studies based on equilibrium thermo- 
dynamics can nevertheless bring valuable information on certain trends in the structural 
evolution of DMSs and with respect to their stability, including the magnetic one. 

The basic information obtained from structural studies based on the effective Ising model, 
e.g., a tendency to clustering, calculated short-range order parameters, etc. can be in turn 
employed in the preparation of better structural model for sampling of magnetic impurities 
in the Heisenberg model as it was demonstrated recently [50!]. We wish to mention that the 
same electronic structure model can be used to calculate the transport properties of DMSs, 
both in the bulk phase and for the multilayer arrangement, which are based on the Kubo- 
Greenwood and Kubo-Landauer approaches as formulated in the framework of the CPA, 
respectively [511 . In this way, the finite-lifetime effects due to various impurities present in 
the sample including native defects and possible co-dopants are included consistently in the 
studies of structural, magnetic, and transport properties of the DMS because they are based 
on the same electronic structure calculated from first principles. 
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